UNCLASSIFIED 


AD NUMBER 


LIMITATION CHANGES 
TO: 


Approved for public release; distribution is 
unlimited. 


FROM: 


Distribution authorized to U.S. Gov't. 


agencies 
and their contractors; 


Critical Technology; MAY 
1968. Other requests shall be referred to Air 
Force Technical Application Center, VELA 
Seismological Center, Washington, DC 20333. 


This document contains export-controlled 
technical data. 


AUTHORITY 
usaf ltr, 28 feb 1972 


THIS PAGE IS UNCLASSIFIED 


4D834 725 


AN ANALYSIS OF A TECHNIQUE FOR THE GENERATION OF 
HIGH RESOLUTION WAVENUMBER SPECTRA 


23 May 1968 


Prepared For 


AIR FORCE TECHNICAL APPLICATIONS CENTER 
Washington, D. C. 


By 
P. R. Lintz 
TELEDYNE, INC. 


Under 
Project VELA UNIFORM 


Sponsored By 
ADVANCED RESEARCH PROJECTS AGENCY 


Nuclear Test Detection Office 
ARPA Order No. 624 


DC 
aren me 
wu UG 
~~ 


"218 


AN ANALYSIS OF A TECHNIQUE FOR THE GENERATION OF 
HIGH RESOLUTION WAVENUMBER SPECTRA 


SEISMIC DATA LABORATORY REPORT NO. 218 


AFTAC Project No.: VELA T/6702 

Project Title: Seismic Data Laboratory 
ARPA Order No.: 624 

ARPA Program Code No.: 8F10 

Name of Contractor: TELEDYNE INDUSTRIES, INC. 
Contract No.: F 33657-68-C-0945 

Date of Contract: 2 March 1968 

Amount of Contract: $ 1,251,000 

Contract Expiration Date: 1 March 1969 

Project Manager: Royal A. Hartenberger 


(703) 836-7647 


P. O. Box 334, Alexandria, Virginia 


AVAILABILITY 


This document is subject to special export controls and each 
transmittal to foreign governments or foreign nationals may 


be made only with prior approval of Chief, AFTAC. len VEL, 
(OC Le333 


a 


YT Ee ort e0 


GRR rar i 


This research was supported by the Advanced Research 
Projecta Agency, Nuclear Test Detection Office, under 
Project VELA-UNIFORM and accomplished under the technical 
direction of the Air Force Technical Applications Center 
under Contract F 33657-68-C-0945 


Neither the Advanced Research Projects Agency nor the 
Air Force Technical Applications Center will be responsible 
for information contained herein which may have been supplied 
by other organizations or contractors, and this document is 
subject to later revision as may be necessary. 


ae 


LIST OF FIGURES 
ABSTRACT 

INTRODUCTION AND THEORY 
PROCEDURE AND RESULTS 
CONCLUSIONS 

REFERENCES 


TABLE OF CONTENTS 


= gy = 


Page No. 


15 
16 


LIST OF FIGURES FIGURE NO. 
A map of the WMO array seismometer grid. 1 


The WMO array response to an infinite-velocity 
plane wave. 2 


The ordinary f-k spectrum of a theoretical 12.0 
km/sec plane wave arriving at WMO from the south; 
the frequency is 0.3125 cps. B 


The high-resolution f-k spectrum of the 12.0 km/sec 
event at 0.3125 cps; S/N = .03. 4 


The high-resolution f-k spectrum for the 12.0 km/sec 
plane waves at frequency of 0.3125 cps. The refer- 
ence seismometer is located at (.94,-.17); S/N = .03. 5 


The high-resolution f-k spectrum for the 12 km/sec 
plane wave at .3125 cps averaged over all of the 
seismometers; S/N = .03. 6 


The ordinary f-k spectrum for the 12 km/sec plane wave 
at .9375 cps. 7 


The high-resolution f-k spectrum for the 12 km/sec 
plane wave at .9375 cps. The reference seismometer 
was at the origin; S/N = .03. 8 


The high-resolution f-k spectrum for the 12 km/sec 
plane wave at .9375 cps. Reference seismometer 
(.94,-.17); S/N = 03. 9 


The high-resolution wavenumber spectrum of the 12 
km/sec plane wave at .9375 cps. The spectrum wae 
averaged over all of the seismometers in the array. 10 


A map of *he LASA subarray Fl seismometer grid. 11 


Ordinary f-k spectrum of LONG SHOr at LASA subarray 
Fl at .625 cps. 12 


High-resolution f-k spectrum of LONG SHOT at LASA, Fl, 
averaged over all seismometers in the array at .625 
cps. C = 2.0. 13 


The ordinary f-k spectrum of LONG SHOT at LASA Fl at 
1.5625 cps. 14 


The high-resolution f-k spectrum of LONG SHOT at LASA 


Fl averaged over all of the seismometers in the array 
at 1.5525 cps S/N= 2.0. 15 


- ii - 


LIST OF FIGURES FIGURE NO. 


The ordinary f-k spectrum of LONG SHOT at LASA Fl 
at 2.50. cps. S/N = 2.0. 16 


The high-resolution f-k spectrum of LONG. SHOT at LASA 
Fl averaged over all of the seismometers in the array 


at 2.50 cps. S/N = 2.0. 17 
Beam-steered power-density spectrum of the LONG SHOT 

event at LASA Fl. 18 
13 sensors from the LASA XE3 subarray. 19 


Array response of the partial XE3 array to an i .*inite 
velocity plane wave. 20 


Ordinary f-k spectrum at .625 cps of two plane waveec, 
one of amplitude 1 at 20 km/sec, North; and one of 
amplitude .2 at 30 km/sec, South. 21 


High-resolution wavenumber-spectrum of the two plane 
waves at .625 cps averaged over all of the sensors in 
the array S/N = 2.0. 22 


High-resolution wavenumber spectrum of the two plane 
waves at .625 cps averaged over all of the sensors in 
the array. S/N = .03. 23 


Ordinary f-k spectrum of the two plane waves at 1.25 cps. 26 


High-resolution f-k spectrum of the two plane waves at 
1.25 cps averaged over all of the seismometers in the 
array. S/N = 2.0. 25 


High-resolution f-k spectrum of the two plane waves at 
1.25 cps averaged over all of the seismometers in the 
array. S/N = .03. 26 


The ordinary f-k spectrum of the two plane waves at 
1.875 cps. 27 


The high-resolution f-k spectrum of the two plane 
waves at 1.875 cps, averaged over all of the seis- 
mometers in the array S/N = 2.0. 28 


The high-resolution wavenumber spectrum of the two 


plane waves at 1.875 cps, averaged over all of the 
sensors in the array. C = .03. 29 


~ iii - 


wey 


ABSTRACT 


In this report, we discuss a method of obtaining a high- 
resolution wavenumber spectrum of time series data obtained 
from a seismic array. We derive the theoretical high-resolution 
wavenumber spectrum of both a single plane-wave event and an 
event made up of the sum of two plane waves of different amplitudes. 
We then process the two theoretical events and analyze the results. 
The plane waves for both events are simulated by impulses, with 
time delays between sensors chosen to simulate move out across 
the array. We also process the LONG SHOT event. 


We conclude that although the high-resolution wavenumber spec- 
trum does indeed give better resolution than the ordinary wave- 
number spectrum, the high-resolution technique does not do a 
much better job of detecting a small signal in the presence of a 
large signal than the ordinary wavenumber spectrum. It appears 
that the problems of detecting a small signal in the presence of 
a large signal might be approached by designing an "approximate" 
k-space whitening filter in the time domain to "pre-whiten" the 
data in k-space, and then +o proceed with ordinary or high- 
resolution frequency-wavenumber analysis. 
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INTRODUCTION AND THEORY 


Frequency -wavenumber spectra 


The ordinary three dimensional frequency-wavenumber 
spectrum (Burg, 1964) of time series data sampled at 
discrete points in space is given by: 

N N 


P(f ,k) = dy ty Sem (w»x.) exp [-i2" k ° (x, ~ =) . 


where an is the mth element of the spectral matrix. The 


fact that P(f,k) is real can be seen by writing it as: 


N N N 
= a + a 
P(f ,k) by $5, (w) + 2 I on [s,,,6w)| 


cos [z= k (x, - x) -@; - “| 3m<j 
where 


8, - on, = wW E (w) - a cw = 27 ko . (x, - x) 
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It is convenient to write the ordinary frequency-wavenumber 
spectrum in vector notations as: 


P(f,k) = AY SA 


where 


+ : é 
A= [exp (-i2m k + X,)>5 exp (-i2n k + x5) *¢¢exp (-i2n k = n> 


is a row vector, the hermitian conjugate of A and where S is 
the N x N spectral matrix: 


The technique of computing high-resolution wave- 
number spectra follows the method used by Haney, (1967). 
The basic method consists of designing a multichannel filter 
whose desired output is an impulse at time t = t, and spatial 
lag x = X5» where t. is usually zero and Xo is the spatial 
position of one of the seismometers. This method, then, 
consists of designing a least-mean-square error "whitening" 
filter to "whiten" the data in f-k Space. The inverse of 
the f-k response of the filter which gives a "white" or 
constant f-k spectrum when applied to the space-time series 
data should give a good estimate of the true f-k spectrum of 
the original data. This method could also be considered as 
equivalent to extending the spatial correlation functions 
(Burg, 1967). 


The-multichannel whitening filter is sensitive to gain 
inequalities. One method of compensation is to normalize 
the spectral matrix (Haney, 1967): 


S.. e 
i - ee 
S s 4% = O55 exp (1055) 


where S55 is the ordinary coherence between channel i and 
channel j and O55 is the phase difference between channel 
i and channel j. Thus the normalized spectral matrix is 


1.0 S45 exp (10,5) *** Gy exp (1054) 


All the information concerning velocity and azimuth is 
included in the coherence and phase, so normalizing the 
spectral matrix at each frequency cannot change the 
apparent velocity and azimuth of an event. Using the 
normalized spectral ratrix, the multichannel filter design 


equations become: 


(er + 6) C= Ty 


where 
id = [o.0, -++2, 0, O, **° o| 
or ‘ 
MU rs Sim 


The high-resolution f-k spectrum designed on the mth sensor 


is given by: 1 
m _4 t+ 1, ns ,\ 
P (fk) -(s H ~H ) 


~ 


“H = (eI + 8)7? ®y = Q Ty 


where ™U is a vector which picks out the mth column of any 
matrix by which it is multiplied. Thus 


Mei” = Mg Mgt = My” 


The spectrum of an impulse at time t=o and position x = x 
and x = o is 


Oo 


E + 2 cos (217 k « x.) | 


- ye 


@ 


This function has a maximum value along lines perpendicular 
to x, and spaced at intervals in k space of k = n/x, . The 
net adil’ of designing a high resolution f-k Speckouit on 
sensors not at the origin but as some % 3 Xo will result 

in a non-physically meaningful stretching oF the f- k spectrum 
along a line perpendicular to Xo° However, calculating the 
f-k spectrum over all of the sensors in the array is useful 
since effects such as seismometer response and ground ccup- 
ling tend to cancel out and the non- Physically-meaningful 
stretching of the f-k spectrum also tends to cancel out for 


symmetric arrays. 


The average spectral matrix of the whitening fiiters 


designed on each of the sensors in turn is 


The quantity 9? is veny easy to calculate, since it is the 
Square of the inverse of some constant times the identity 
matrix plus the smoothed Spectral matrix. The constant (c) 
can be considered to be the signal-to-noise ratio of the 
random signal spectral matrix I. 


Response of the high-resolution technique to a single 
etn technique to a single 
plane wave input 
—_———$______ 


In the special case of an impulsive plane wave with unity 
conerence between sensors, the spectral matrix may be factored 
as follows: 


+ 
where u 


but 


+ 
and u 
so that 


F + 
Since A 


co 

cS 

c 
D> 


is the row vector 


[exp (-iwt,), exp (-iwt,), *** exp (-iut)] 


-1 
(cI + uu*) = 2 aC. )=Q 
ctu u 
= N 
2 _ 1 uu N 
i ae [ - cen © c+ 7? | 
IA=zN 


E At + A 
+ 2 am 1 ~ UuU wn N 
4 Q A = a - S07 (9 = 3D] 


= P(f,k) = ordinary f-k spectrum 


(c + n)71 may be expanded for c<N as foliows: 


HR (f,k) = 
t 


ee 


N N N 
P(f,k) = J Sz.(w) + 2 J L |8z_¢w)] cos [2x (k-k 9 )* (x. =x | 
jon 1 j=2 m=. JM md Mac Mi 
m<j 
at k = ko 
P N _ N(N-1)] _ 12 
P(f,k) * Fag = NS -2C, SMe 2 racy] 
Thus 


oN 
HR (f,k)) = 5 
Cc 
n-w h -€2)] 


hence the maximum value of ®HR(f,k) is independent of c up 


to terms of the order ; 


= NA 


when 


P(f,k) + 0 


OHR(£,k) + 07 


Response of the high resolution technique to an input 
consisting of two plane waves input 


In the case of where the event consists of two plane 
waves of unity coherence, one wave of amplitude unity, the 
other of amplitude b, the ordinary f-k spectrum is: 


+ + + 
P(f,k) =A (uy + bu,) (uy + bu,) A 


wae: eee ae oe: Nae, 


Let 


then 


where 


aa 


P(f,k) = P,(£,k) + b° Pi(£,k) + 2b A* GCw) A 


G. 
1 


5) = 


cos [enece; - 7) 


where ti is the arrival time of piane wave number 1 at sensor 
i and ts is the arrival time of plane wave number 2 at sensor 


Acs 
For the high-resolution frequency-wavenumber spectrun, 

1 (u, + bu,) (uy + bu, ) 
Praem) ed f SES 
él + Cui. * ba.) =Q:— 

1 2 c c + (u, + bu, ) (uy + bu,) 
Ca; ® bus) Gi, + buyd= C2 + b4) we BOT de out aw ) 
1 2 1 2 12 2 1 
N ing 
b(u,tu Uy u,) = 2b Qh cos [2ne (ts - tp) = 2b Tr [eww] 
+ + 
ee: : (u; + buy) (u, + buy) 
Cc ec + Rw) 


PE 


OH 


where 


R(w) = (1 + b?)N + 2b 1 [ow] 


ae = (uy + bu, ) (uj + bu, ) _ _RWw) 
U8 ee ee S + RUD 


i> 
vie 


to PCE »k) _ Rw) 
c + Rw) c+ R (Ww) 


and where 


P(f,k) is the ordinary f-k spectrum 


P(f,k) = Py (fk) + bP, (fk) + 2b A* Glw) A 


OHR(f,k)= ——N___ 


b> 
2 
b> 


-1 
Expanding E + Rw | and simplifying, we have 


OHR(f,k) = 


c oN 
"HPD Gta) Gia) 


The expression for the high-resolution wavenumber spectrum 
for two plane waves input is seen to be quite similar to the 
expression for the spectrum of a single plane wave, with 

the exception that F(f,k) has cross terms, and N is replaced 
by R(w) in the expansion (c + n)7!, 
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PROCEDURE AND RESULTS 


Procedure 


Synthetic spike seismograms were generated for the WMO array 
and for the LASA XE3 array. The WMO seismogram corresponded to a 
12 km/sec impulse arriving from due south of the array. The XE3 
seismogram consisted of the sum of two impulses, one from due north 
with an apparent horizontal velocity of 20 km/sec and an amplitude 
of unity, and the second impulse from due south with an apparent 
horizontal velocity of 30 km/sec and an amplitude of 0.2. The 
LONG SHOT event at LASA subarray Fl was also processed. The 
spectral matrix of each input was calculated and smoothed using 
the fast Fourier transform program of Claerbout, et al. (McCowan, 
1966). Both ordinary and high-resolution spectra were calculated 
for each event. High-resolution spectra calculated on individual 
sensors and averaged over all sensors (as discussed in Section A) 
were generated for the single plane wave at WMO, while only the 
"average" high-resolution spectra were calculated for the two-wave 
event at LASA XE3. 


The ordinary and the high-resolution f-k Spectra of a single ' 
plane wave input 


Figures 3 and 7 are the ordinary f-k spectrum of the 12 km/sec 
impulse at frequencies -3125 cps and .9375 cps. Figures 4 and 8 
are the high- resolution wavenumber spectrum (designed on the center 
seismometer) for the same event at the same frequencies. Comparing 
Figure 4 with Figure 3, we see that the high-resolution wavenumber 


spectrum at .325 cps indeed shows an improvement in resolution over 


_ the ordinary spectrum (Figure 3). Figure 4 shows a signal arriving 


from the south with a velocity of slightly greater than 10 km/sec, 
whereas Figure 3 (the ordinary f-k spectrum) shows a Signal coming 
from the general direction of the south but with much uncertainty 

as to the exact velocity and azimuth. At 0.9375 cps, the ordinary 
f-k spectrum (Figure 7) shows a plane wave from a southerly direc- 
tion, but the amount of uncertainty as to velocity and azimuth i 


- ll - 


large compared to tha: of the high-resolution spectrum (Figure 8). 
Figures 5 and 9 are the high-resolution wavenumber spectrum (de- 
signed on the sensor located at (.94, - .17) at frequencies .3125 
cps and .9375 cps. Comparing Figure 5 with Figure 4, we see that 
the resolution of the high-resolution f-k spectrum (Figure 5) de- 
Signed on the seismometer at (.94, - .17) is approximately equal 
to the resolution of the high-resolution f-k Spectrum (Figure 4) 
designed on the center (0,0) seismometer; however, the spectrum 

in Figure 5 is "stretched" along the line perpendicular to the 
direction of the line connecting the origin to the position 

(.94, - .17). In other words, the direction of the stretch is 

90° + tan71 (- .17/.94) = 90° - 16°, or 9 = 74°, The resolution 
of Figure 9, the high-resolution spectrum at .9375 cps designed 

on the off-center seismometer is approximately equal to that of 
Figure 8, but again the Spectrum is stretched in the direction of 
74°, Thus the net result of designing a high-resolution f-k spec- 
trum on a seismometer located at x= xX, is indeed to introduce a 
cos (2mk+x, ) term, which causes a stretching of the spectrum along 
a direction perpendicular to the vector x = xX). Figures 6 and 8 are 
the high-resolution f-k spectra averaged over all of the seismo- 
meters in the array, for the 12 km/sec event, at frequencies .3125 
cps and .9375 cps. Comparing Figure 6 with Figures 4 and 5, we 
see that in the high-resolution spectrum which has been averaged 
over all of the seismometers, the non-physically meaningful 
stretching of the spectrum due to designing the high-resolution 
Spectrum on an off-center seismometer has been cancelled out. 

The spectrum of Pigure 5 seems to be not quite as peaked as that 
of Figure 4, although the difference is insignificant. The spec- 
trum of Figure 4 is also tilted to the left due to the asymmetry of 
the WMO array, whereas the spectrum of Figure 5 seems to be Slightly 
"stretched" in a northwardly direction. Figure 10 is the high- 
resolution wavenumber Spectrum averaged over all seismometers in 
the array, at .9375 cps. The averaged spectrum has slightly less 
resolution than the spectrum that was designed on the ‘center seis- 
mometer (Figure 8); however, the "stretching" of the spectrum 
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evident in the spectrum designed on the off-center seismometer 
(Figure 9) has essentially been cancelled out. The averaged spec- 
trum should give better results when effects such as different 
seismometer responses and ground coupling for each sensor are con- 
sidered. 


Figures 12, 14, and 16 are the ordinary f-k spectrum of LONG 
SHOT recorded at LASA Fl at frequencies 0.625, 1.5625, and 2.50 cps. 
Figures 13, 15, and 17 are the "averaged" high-resolution f-k 
spectra at the same frequencies. The signal-to-noise ratio (c) was 
2.0. At 0.625 eps, both the averaged high-resolution (Figure 13) 
and the ordinary spectrum (Figure 2) indicate the event is coming 
from the correct direction, but the velocity is somewhat smaller 
than expected (11.6 km/sec, vs 14.1 km/sec, from the J-B travel- 
time tables). At 1.5625 cps the high-resolution spectrum (Figure 
15) shows a velocity of approximately 18 km/sec, whereas the 
ordinary f-k spectrum, (Figure 14) shows a velocity closer to the 
velocity computed from the travel-time tables, although the reso- 
lution is such as to leave a large uncertainty as to the exact ve- 
locity and azimuth. At 2.50 cps both the high-resolution (Figure 
18) and the ordinary f-k spectrum (Figure 17) agree quite well as 
to the velocity and azimuth of LONG SHOT at LASA Fl. The velocity 
of the event is quite close to the predicted 14.1 km/sec velocity. 
The high-resolution spectrum offers a little more resolution and 
hence less uncertainty than the ordinary f-k spectrum. Figure 18 
is a beam-steered (time-shift and sum) power density spectrum of 
LONG SHOT. Nine seismometers in-line in the approximate direction 
of the event were time-shifted and summed. The power spectrum of 
the sum was then computed and displayed in Figure 18. This demon- 
strates that there is indeed enough power at 2.5 cps to give an 
accurate f-k spectrum at that frequency. 


The ordinary and the high-resolution f-k spectra of two plane 
t 


Waves inpu 


Figure 20 is the array response of thirteen of the seismometers 
at LASA XE3. Figures 21, 24, and 27 are the ordinary f-k spectrum 
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of a plane wave of velocity 20 km/sec and amplitude 1.0 plus a 
plane wave of velocity 30 km/sec and amplitude 0.2, computed at 
0.625, 1.25, and 1.875 ops. Figures 22, 25, and 28 are the aver- 
aged high-resolution f-k spectrum (S/N = 2.0) of the same event 
calculated at the same frequencies. Figures 23, 26, and 29 are 
the high-resolution f-k spectrum (S/N = 0.03) of the same event, 
calculated at the same frequencies. At .625 cps the 30 km/sec 
event would not be detected by either the ordinary f-k spectrum 
(Figure 21) or by the average (S/N = 2.6) high-resolution f-k 
spectrum (Figure 22). The extreme elongation of the average 

(S/N = .03) high-resolution f-k spectrum (Figure 23) might make 
One suspicious about the presence of a second event, but the con- 
clusion would depend on the subjective discrimination of the ana- 
2yst. At 1.25 cps, neither the ordinary f-k spectrum (Figure 24) 
nor the averaged high-resoluiion (S/N = 2.0) £-k spectrum offer 
much possibility of detecting the 30 km/sec plane wave in the 
presence of the larger 20 km/sec event. The averaged high- 
resolution (S/N = 0.3) f-k spectrum offers a distinct possi- 
bility of detection, although the resolution of the two peaks is 
far from optimum. At 1.975 cps, the 30 km/sec event would prob- 
ably be mistaken for a side lobe of the ordinary f-k spectrum 
(Figure 27). The averaged high-resolution (S/N = 2.0) f-k spec- 
trum (Figure 28) does not give any indication of the 30 km/sec 
event. The averaged high-resolution (S/N = .03) f-k spectrum 
gives a definite indication of the smaller event, but there is 

a spurious 6 db peak at 10 km/sec from the southwest. The com- 
plicated spectra of an event made up of a sum of two plane waves 
is due to the cross-terms 4*G(w)d4 in the frequency-wavenumber 
Spectrum which give rise to spurious peaks. 
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CONCLUSIONS 


From this report, we conclude that although the high-resolution 
frequency-wavenumber spectrum gives much better resolution for 
single events, the high-resolution technique does not improve sig- 
nificantly the Capability of a seismic array to detect multiple 
time-overlapping events from different azimuths. A possible method 
of detecting a small event in the presence of a large event is, we 
Suggest, to design an approximate pre-whitening filter in the time 
domain to approximately "whiten" the f-k spectrum of the data, 
convolve the filter with the time series data and then proceed with 
either ordinary or high-resolution f-k analysis. The fact that 
the pre-whitening filter need be only an approximate whitening 
filter would leave quite a bit of lattitude in its design. 
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Figure 2. The WMO array response to an infinite-velocity 
plane wave. 
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Figure 3. The ordinary f-k spectrum of a theoretical 
12.0 km/sec plane wave arriving at WMO from 
the south; the frequency is 0.3125 cps. 
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Figure 4. The high-resolution f-k spectrum of the 12.0 
km/sec event at 0.3125 cps; S/N = .03. 
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Figure 5. The high-resolution f-k spectrum for the 12.0 
km/sec plane waves at frequency at 0.3125 cps. 
The reference seismometer is located at 
(.94,-.17); S/N = .03. 
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Figure 6. The high-resolution f-k spectrum for the 12 
km/sec plane wave at .3125 cps averaged over 
all of the seismometers; S/N = .03. 
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Figure 7. The ordinary f-k spectrum for the 12 km/sec 
plane wave at .9375 cps. 
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Figure 8. The High-resolution f-k spectrum .or the 12 
km/sec plane wave at .9375 cps. The reference 
seismometer was at the origin; S/N = .03. 
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Figure 9. The high-resolution f£-k spectrum for the 12 
km/sec plane wave at .9375 cps. Reference 
seismometer (.94,-.17); S/N = .03. 
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Figure 10. The high-resolution wavenumber spectrum of the 
12 km/sec plane wave at .9375 cps. The spectrum 
was averaged over all of the seismometers in 
the array. 
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Figure 11. A map of the LASA subarray Fl seismometer grid. 
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Figure 12. Ordinary f-k spectrum of LONG SHOS at LASA sub- 
array Fl at .625 cps. 
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Figure 13. High-resolution f-k Spectrum of LONG SHOT at 
LASA, Fl, averaged over all seismometers in 
the array at .625 cps. C= 2.0. 
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Figure 14. The ordinary f-k spectrum of LONG SHOT at LASA 
Fl at 1.5625 cps. 
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Figure 15. The high-resolution f-k spectrum of LONG SHOT 
at LASA Fl averaged over all of the seismometers 
in the array at 1.5625 cps C = 2.0. 
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: _ Figure 16. The ordinary f£-k spectrum of LONG SHOT at LASA 
Fl at 2.50 cps. C = 2.0. 


a 


KK WAVG NUMUER WG Mak, Value © FREQUENCY & 2.50000 


o2,359E201 =2.00%"n1 oly? ff 2 & S.00F+o2 1-o0F 
K E.BORelecne 46666606666660066606 
vy | 666606666666666066 
1 140666646606666566666: 
C] 1 I 66666666466666666 
A 1 I 66666666000666 
v 1 1 $666660666660066006 
8 Fe BOG+9l tense enn ten, ncneettacecnon ns $set ete Suet sen heneens con c00h 6066066666666 7600 -nne te ne teeeectaue ~=966t 
! 1 ! ‘ i] i] 4666666666660 i] i] 6666! 
® ' 3 I 06660606666 1 om | 506666) 
v ! 1 | 666066066 i] 1466666661 
" ! 1 1 666066666686 1 6666665666661 
s ! i) i 666 0066666606066606666 66666666666666661 
© 1,808+01 +--+" -+-~ enon Geer ngecengcoce=08 ©0416 60606666060066666666 66066666606060666664 
® ! 1 1 606 66066066600666656660666666666666666i 
Q Boy 1 1$60666666666666666606609066666666666666) 
J 1 1 186666666666666666566666%6666666666666351 
! 1 i] 1 6666666666660666666666%6666666666666.\61 
! 1 1 1 666666606666666666666%666666 6666666641 
1. 008-G1¢+-*+--F0000080Q00U YoUNQO? UBUD Sed Pwee oe Cog conces 2006009 $66666666666666666669666664666666666+ 
! 1 1 1 66666666666666666 66°66 6666666666666) 
‘- - -t. 3 1 - — - $666666666666666606666666666666661 
' 1 1 i] * $6666660666666606666666666666661 
! 1 1 1 J 664e-1 1 66666666! 
' 1 1 1 J i] 1 1 
PT ee eM eB 00 he ccee cote ftec ces cee bow cc ccces tool ween etone wnt eeheconvnte ct 
! 1 1 1 1 1 i} 1 1 
+ i i] + . 1 + 1 t 1 1, I 
' ! ! ! 1 1 1 1 1 1 1 
! ' ' i] 1 1 1 1 1 1 1 
! ' 1 1 1 1 I 1 1 1 1 
08 COB O91 Oe we cer cree tegen cece 0% $0 208 08 8G $ S059 00h 9 08SR cel eh Scene Fe0e bowen sce $2 wbasuned $nceFenceebocencnec ct 
1 1 1 1 1 1 I I 1 1 
t at -I- + ao a | See | - 1 1 1 
| 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 d 1 i] 1 
1 1 1 1 et. 1 
05 298192 6 en we cece $000.02 099% $ eo Fe Soe 0 8 $99 ES eg Pooh ncesccecgccconyCodbucceccces } '$66666666666* en mvetoenncncect 
1 1 1 1 1 1 666666666666066066 = 1 I 
§—--—-.-. 4. + 4. 7k - §- - we deo 6066 -4 
MEETEEVOCS 1 1 1 1 1 1 
066 666666660686 1 1 1 1 1 6 i] 
aseesescausnsasedy | 1 1 1 1 66 1 
4.008 -0166066060060 0B06 oc ccccrccgooes cotoctBoune cevetecccvccccto enh bh = 
+ 066660606660 6008 1 1 ] 1 
666066666600 6006 +. 1. : 666666666650666666066666663 
66666000666066666600666' 1 3 1 6666666666660! 
666666066660666060066666 1 4 1 6 6 
6660666066606666000004606 1 1 1 6 6661 
91.308 =0366606060666 06660606066 606 bho teens cn tocpaccopccgegooccce$h bbb 6 1666664664 
666660660660606666606660666' 1 1 i] 6 6666661 
6666606066606066006060606666a. 3. 1. . -b. --+66 6666661 
1 | 666666666661 
666066606660066066606666666606) 1 66666666666) 


6660606066606 6606660666066666)) 
66606666660066006666666666666' 


1 
1 1 1 
oP oper ehos soccegerergos: 


©2008 -9166660600666006606060646066666 ee 
6660666666006650666066666606061 1 1 61 
00660006666060006000666060006 1 L- iE rt 
€6406066666066606000666606666 | 1 1 61 
06666060b0b0bK60bb00EC008000 1 1 1 ‘* 61 
06606666666066660600666600668 | 1 1 466 61 


92 Bob egl tone ccwccesrcccce ne *qamgpapee= $oenes cb ovdveccocetoteccsccccetocce Ore ccccocetoce oe 


Figure 17. The high-resolution *-k spectrum of LONG SHOT 
at LASA Fl averaged over all of the seismometers 
in the array at 2.50 cps. C = 2.0. 
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Figure 18. Beam-steered power-density spectrum of the LONG 
SHOT event at LASA Fl. 
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13 sensors from the LASA XE3 subarray. 
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Figure 20. Array response of the partial XE3 array to an 
infinite velocity plane wave. 
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Figure 21. Ordinary f-k spectrum at .625 cps of two plane 
waves, one of amplitude 1 at 20 km/sec, North; 
and one of amplitude .2 at 30 km/sec, South. 
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Figure 22. High-resolution wavenumber-spectrum of the two 
plane waves at .625 cps averaged over all of 
the sensors in the array S/N = 2.0. 
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Figure 23. High-resolution wavenumber spectrum of the two 
plane waves at .625 cps averaged over all of 
the sensors in the array. S/N = .03. 
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Figure 24. 


Ordinary f-k spectrum of the two plane waves 
at 1.25 cps. 
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Figure 25. High-resolution f-k spectrum of the two plane 
waves at 1.25 cps averaged over all of the 
seismometers in the array. S/N = 2.0. 
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Figure 26. High-resolution f-k spectrum of ‘the two plane 
-waves at 1.25 cps averaged over all of the 
seismometers in the array. S/N = .03. 
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Figure 27. The ordinary f-k spectrum of the two plane 
waves at 1.875 cps. 
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Figure 28. The high-resolution f-k spectrum of the two 
plane waves at 1.875 cps, averaged over all 
of the seismometers in the array S/N = 2.0. 
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Figure 29. The high-resolution wavenumber spectrum of the 
two plane waves at 1.875 cps, averaged over 
all of the senso~s in the array. C= .03. 
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